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Abstract:- alternating direction explicit and alternating direction implicit methods (ADE and ADI) were used 
to solve Schnackenberg model, we were found that alternating direction implicit method is much more accurate 
and faster than alternating direction explicit in this kind of models. 
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I. INTRODUCTION 

Reaction-diffusion (RD) systems arise frequently in the study of chemical and biological phenomena 
and are naturally modeled by parabolic partial differential equations (PDEs). The dynamics of RD systems has 
been the subject of intense research activity over the past decades. The reason is that RD system exhibit very 
rich dynamic behavior including periodic and quasi -periodic solutions [6]. 

Various orders are self-organized far from the chemical equilibrium. The theoretical procedures and 
notions to describe the dynamics of patterns formation have been developed for the last three decades [4]. 
Attempts have also been made to understand morphological orders in biology [5], Clarification of the 
mechanisms of the formation of orders and the relationship among them has been one of the fundamental 
problems in non-equilibrium statistical physics [3]. 

I.I. MATHEMATICAL MODEL 

A general class of nonlinear-diffusion system is in the form 
du > 

— = d L Au 4- a L u 4- £) L i? -I- /(iia") + SiGc3 
at 

dv 

— = d 2 Au +- c : u -I- b 2 v -f(u..v) 4- 5 : Of] t 
at 

With homogenous Dirchlet or Neumann boundary condition on a bounded domain Q. , n<3, with 
locally Lipschitz continuous boundary. It is well known that reaction and diffusion of chemical or biochemical 
species can produce a variety of spatial patterns. This class of reaction diffusion systems includes some 
significant pattern formation equations arising from the modeling of kinetics of chemical or biochemical 
reactions and from the biological pattern formation theory. 

Schnackenberg model: 

a L = — k, b L = q 2 = b z = 0, / = a 2 v,g l = n,g% = h, 
where k, a and b are positive constants. 

Then one obtains the following system of two nonlinearly coupled reaction-diffusion equations (the 
Schnakenberg model), 

— = d t Au - ku + u z v + a, (t, x) G (0, oo] x £2 
at 

dv 

— = d-Av -u 2 v + b, (t, x) E (0, «s] x SI 
at 

u(t,x) = vlt.x) = 0, t>0, xsd£l 
u(0, x) - u (x}, v(0,x) - v Of] , xeSl 

Where d L ,d 2 , a, k and b are positive constants [8]. Various finite difference algorithms or schemes 
have been presented for the solution of hyperbolic-parabolic problem or its simpler derivatives, such as the 
classical diffusion equation. It is well-known that many of these schemes are partially unsatisfactory due to the 
formation of oscillations and numerical diffusion within the solution [1,7]. 

Solution by the finite difference method, although more general, will involve stability and convergence 
problems, may require special handling of boundary conditions, and may require large computer storage and 
execution time. The problem of numerical dispersion for finite difference solutions is also difficult to overcome 

m. 
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II. MATERIALS AND METHODS 



1 1.1. Derivation of alternating direction explicit (ADE) for schnackenberg model: 

The two dimensional Schnakenberg model is given by 



4- ji'v 4- a.. 



01- 
~3t 



dx 2 dy 1 



(t.x) t (Q.-i) x Q 

(t, x) G (0, oo) X £1 



(2) 



We consider a square region 0<x<l, 0<y<l and u, v are known at all points within and on the 
boundary of the square region. We draw lines parallel to x, y, f-axis as x=\h, y=]k, and t=nz, i,j=0, 1, 2,...,M 
and n=0, \,2,...,N, where h=Sx, k= dy, z= St. 

The explicit finite difference approximation (ADE) to shnackenberg model in two-dimensions are given by: 

Y t _ d i I £5 } + d ^ \ p J ~ k "*» + 

+ a 



5t = ^ I ^ J + da I fc*" } ~ lff W ^ + * 

Multiplying both equations by and set fe = k, then we have a rectangular region and replacing ~rf- = r L and 

—7- = r- ,then we get 

^L.j'.n+i — ^ij'.R = 

*1 (Pi-UjJi - 2U i,j,n + Vi + L.j,n + U i.;-l,Ti - 2U i r j,n + ^1,;+ 1») ~ M£ ^ij.n +" 
and 

Then simplifying the system to obtain 

= (1 - 4r L - kSt)U i:i i , +■ r L (^_ Un +■ t/ £+u _ B 4- t/ u _,, n + (J u+1 ,J + fft (l/yjVy. „ + oft 

n + l = (1 - +" ft-U.n + <W + Kj-u, + C iJ+lj J - ^ Oijj\ M + 5t b 

This is the alternating direction explicit formula for the Schnakenberg model 



II.2. Derivation of alternating direction implicit (ADI) for schnackenberg model: 

In the ADI approach, the finite difference equations are written in terms of quantities at two x levels. 
However, two different finite difference approximations are used alternately, one to advance the calculations 
from the plane n to a plane n+1, and the second to advance the calculations from («+7)-plane to the (n+2)- 

plane by replacing — - and — - by implicit finite difference approximation [5]. we get 

st = M — ~ * ) + d A ^ )- kU i,.« 

J t = d 2 { # J + d 2 { ^ J - Wj,iJ + 

we set h = fc, then we have a rectangular region and multiplying the equations by St and replacing -rf^ = i\ 
and -r^ - = f 2 , then we get 

+ ^ {U iiipt )\^ + £ta 

Vi,j,n+X ~ = ^fa-l.,;.,P. + L - ^L.J..n. + l + ^1.^1 + ~ 2^.,;.,n + Kj+lj.) ~ (fi,;,)^,;* + fft * 
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r 2^i-l..j,n+L + 2*2^L,;„n+L *2^i 1,11+1 

= r 2 Vi + & - 2^)1^ + r 2 If J+Wl - St (l/^J 3 ^ + -fft Ji 
+ Cl - Zltyj- n+ i - ^i + i.j. n+ i = (1 - 2l^ ti ,„ + 11^^ + flta + St (t/ u ..j\ M 

or 

— l^i-lj.n+l + (l — 2r 1 )[/[ : j p _ + 1 — l^i + lj'.n + l 

= (1 -2r t )U iln +■ r L I7 u+Ln + fft U .n)Vj.n + ft (o - fcl/yj 

and 

— ^Z^-l.j.n+l + (l — 2r a )F i j j1l+L — +l, J f+l s n+l 

= ^.j-i* +C1- 2r 2 )lf.^ + r 3 lf J+w - ff£ (ffyjX,-,, + fit bt 
Advance the solution from the (n+l) th plane to (n+2) th plane by replacing —7 and — — with explicit finite 

difference approximation at (n+l) th plane then j-t and y-7 by an implicit finite approximation at the (n+2) th 
plane 

n + : — "i,/ n+1 /^i-l j'n+l — — j' n + 1 + ^L + l.,/,n + l\ f^i.j-lji+z — ^"e,;,h+2 + "ij'+lji+aX 

a = ^ 1 ^ J + d 4 p ) 

and 

ffi = d 2 I £i ) + d * I. J - W.j.nJ *U" + fi 

Multiplying the equations by St and replacing r t = -tt~ and r 3 = when h = k, then we have 

^i,;,n+2 ~ ^i,/,Ti+l = 1 (^i-l^n + l — ^ijiti+l + ^l + 1,/,ti+i) + t i (^i,;-Ui+2 — 21/(^+2 + ^1,^ + 1,11+2} 

- kSt U UJ] + it Oij,j\;, n + 5* D 

and 

^LJ.Tl+Z _ ^J,fl+1 = 'zt^L -l.j>. + l _ 2F(j ifl+ ]_ + I^+^jj+l) + 1*2(^1,/- 1,71+2 _ 2^j ifl+2 + ^i,; + 1,11+2) 

and this implies that 

— 1^LJ-l,Tl+2 + (l + 2f£ ^jj^a ~ r L^i,j + Un+2 

= ^i-lj.n + (1 - 2?iM.j. + i n +" ^i + lj+.n + l + « («~ fcU i,j,J + fa.>,)\,, 

and 

~ r 2^i,j-ua+2 + (l — 2? a )l'' i j yrl+2 — r z V{ j +ljfl+ 2 

= T"2^i -ij>- +1 + (1 - 2r : 3 ^ + V i+Un+1 - fit (l/ij >B ) 2 ^ J>T1 + fit i 
Then we get two systems 

— T l^LJ-l,Tl+2 + (l + 27i.]l/y jT1+ 2 — r l^LJ + l,Tl+2 

= ii£/ £ _ u . n + (1 - 2^ j+1 n +■ r L y i+1 J+ _ B + 1 +5t (a - kV Uj J + (u irj j\ jjL 
~ r 2^i,j-ua+2 + (l — 2r 2 )F i j j11+ 2 — la^.j+^H+a 

= ^i-lj>.+l + (1 - 2 ^ V ij.n + l + Iit5+lj**l - ft faj,n)Vj.n + ft & 

and 

— I^L-lj.R + l + (l — 2rj_]lf £ j n + 1 — 1^L + l,;„TWl 

= (1 -2^U; Sn +r t U iJ+Un + St (U un fv ijn + St(a- kU^J 
~ r 2^i-i..j,n+L + tl — 2r 2 )F i j j . rl+ ]_ — ia^+i.j+iji+1 

= ^L,j-m + d -2^)^^ +7,Vi iMl! -if fejj^+fiii 

The last two systems represent Alternating Direction implicit method under the conditions 
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u l,q,n+\ = u l,q+l,n+l 



11 M ,q,n + \ ~ U M ,q+l,n+\ 



and 



v 



\,q,n+\ 





= 0. 



M ,q,n+\ v M ,q+\,n+\ 

The tridiagonal matrices for the system in the level n advanced to the level n+1, for both u and v can be 
formulated as follows AU=B. 

u 2,q,n+l 
u 3,q,n+l 
u 4,q,n+l 
u 5,q+l 



+ 2rj) 


-1 




















"1 


(1 + 2^) 


-1 























-1 


(l + 2q) 


-1 























"1 


(l + 2r x ) 


-1 


















"1 






































































-1 


(1 + 2/j) 


-1 























"1 


(l + 2rj) 



">i-3,g+l 
u M-l,q,n+\ 



■ + r 2 ( u 2,q+\,n + u 2,q-\,n ) + (1 - 2r 2 - z)w 2>9> „ + 2,q,nV 2 , q ,„ 

2 

• + r 2( M 3 lI? +Ui + "3,<z-l,n) + C 1 - 2r 2 - z)"3 ><7 ,« + z " 3, ? ,nv 3) ^„ 

2 

' + r 2 ("4,4+l,« + "4,<7-l,« ) + (! - 2r 2 - z)«4, 9 ,« + z" 4^,nv 4i?; „ 



az + r- 



'2 + u M-\,q-\,n ) + (1 - 2r 2 - z)« M -l, 9 ,n + z " M-l, ? ,nV M _ 1>? n 



d + 2<i) 


-'1 




















"I 


(l + 2/l) 


-1 

























(1 + 2/1) 


-1 























"1 


(1 + 2/1) 


-1 


















"'I 






































































-1 


(l + 2/l) 


-\ 























-1 


(l + 2r 1 )_ 



v 2,4,n+l 
v 3,q,n+\ 
v 4,q,n+\ 
v 5,q,n+\ 



v M-3,q,n+l 
v M-2,q,n+l 
v M-l,q,n+l 
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>2 ( v 2, 9 +l,n + x Xq-U ) + (! " %r % ~ zK *> v 2,q,n - *» U^2,q,n 

2 

"zHq+U + v \q-\,n) + C 1 - 2r 2 ~ ^ v 3,q,n - zu Xq,n v 3,q,n 

2 

1.(%q+l,n + V 4 >9 -1,«) + ( ! - 2r 2 ~ zK H,q,n ~ zu ^q,» v A,q,n 



l( v M-l,q+l,n + v M-\,q-\,n) + (l ~ 2r 2 ~ zK ~> v M-l,q,n ~ zu ^-Kq,n v M-l,q,n 

And the tridiagonal matrices for the system in level n+1 advanced to level n+2 for both u and v are 
given by AU=B. 



o 



(l + 2r 2 ) -r 2 

-r 2 (l + 2r 2 ) -12 

-r x (1+2)!) -ij 

-r 2 (1 + 2r 2 ) 

-r 2 











-r 2 (1 + 2*2) ~ r 2 
-r 2 (1 + 2r 2 ) 



"p,2,n+2 
u p,3,n+2 
M p,4,n+2 
"p,5,«+2 



"p,M-3,«+2 
u p,M-2,n+2 
u p,M-l,n+2 



2 



/X + r l(>+l,2,n+l +"p-l,2,«+l) + < 1 - 2r l - *")«p,2,«+l " 2 "p,2,n + P,%nv p Xn 

2 

^ + r l (";>+l,3,n+l + "p-l,3,n+l ' + " 2 1 _ z >"/),3,n+l " z "p,3,n + z " PA" v p,3,n 

2 

^ + l("p+l,4,n+l + " /?-l,4,n+l ) + (1 - 2r l -z)"p,4,n+l - z 'V,4,n + z " PA nV pA,n 

2 

^ + 1 ("p+l,5,n+l + »p-l,5,n+l ' + " 2, 'l _ : >"p,5,n+l " z "p,5,« + » pM^p,5,n 



fX + r l("p+l,M-\,n+\ + "p-l,M -l,n+0 + ( - l - 2r l -^ u p,M-l,n+l-~ u pM-l,2,n + pA»* 'pM-l,«+l 

Also for AV=B the tridiagonal is in the form 



l + 2m 2 ) 


~ m 2 




















-m 2 


(l + 2m 2 ) 


~ m 2 























-m 2 


(l + 2m 2 ) 


~ m 2 























_m 2 


(l + 2m 2 ) 


~ m 2 


















-my 






































































"«2 


(l + 2ra 2 ) 


-m 2 























-m 2 


(l+2m 2 ) 



v p,2,n+2 
v p,3,n+2 
v pA,n+2 
v p,5,n+2 



v p,M~3,n+2 
v p,M-2,n+2 
v p,M-\,n+2 
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2 

M l(VU»+l +v P-lA»+l ) + (1 " 2 ' n r^VM +zu p,2,n "» Pt* v pXn 

2 

m l < v p+l,3,n+l + v p-l,3,«+l ) + (! - 2m l - z*O l >,3,«+l + z« p>3j „ - z« p,3, n v p 3 n 

2 

m l( v p+lA,n+l +v P -l,4,n+0 + ( l - 2m l-z K >pA,n+l + z u pA,n- zu pA* v p,4ji 



m \ ("p+Uf -U+l + v p-l,M-l,n+l ) + - 2m l - zK ) v p,M-l,n+l + zu p,M-ln ~ zu PM-lnv p M _y n 



III. Application (numerical example) 

Example: We solved the following example numerically to illustrate the efficiency of the presented methods, 



suppose we have the system 




We the initial conditions 

uix,y, 0),v(x,y, 0)) = (exp(-j: - y) , expU + y)) 

y, t), v(0, y, t» = (exp (-t/2 - y) ,exp(t/2 + y)) 
«Cx, 0.. t), v(x, 0, t]) = (exp(-t/2 - x) .. exp (t/2 -f x)) 



[u,!?'}') = (exp (— — x — yj , exp (~ + * + yj) is the exact solution of the problem. 
Where a=b=0.1. 

Then the results in more details are shown in following table and figure: 



ADE-ADI 



3 




0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



0.3 



space 
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Fig. 1: Comparison of ADE, ADI with Exact solution (u). 



Table 1: Comparison of ADE and ADI with Exact solution (u). 



X 


ADE 


ADI 


Exact 


0.1 


2.69265 


2.6912344 


2.6912344 


0.15 


2.423622 


2.27384012 


2.20573840 


0.2 


2.181472 


2.01678803 


1.8971263 


0.25 


1.963516 


1.8279602 


1.6823796 


0.3 


1.767336 


1.64066 


1.52614066 


0.35 


1.590757 


1.4899208 


1.401599208 


0.4 


1.431821 


1.31680 


1.254730168 


0.45 


l. 288764 


1.18920286 


1. 14620286 


0.5 


1.160001 


1.0435533 


1.011435533 


0.55 


l. 044102 


0.939289 


0.879392894 


0.6 


0.939783 


0.8454428 


0.78454428 


0.65 


0.845887 


0.7609726 


0.697609726 


0.7 


0.761373 


0.6849420 


0.635849420 


0.75 


0.685302 


0.6165078 


0.568960165 


0.8 


0.616832 


0.5549110 


0.52549110 


0.85 


0.5552032 


0.4994686 


0.44994686 


0.9 


0.4997315 


0.4495655 


0.40495655 


0.95 


0.4498021 


0.4046483 


0.374046483 


1 


0.4048613 


0.403642189 


0.363642189 



IV. CONCLUSION 

We saw that alternating direction implicit is more accurate than alternating direction explicit method 
for solving Schnakenberg model. 
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